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Abstract 

The aim of the paper is to propose geometric descriptions of multipartite entangled states using algebraic 
geometry. In the context of this paper, geometric means each stratum of the Hilbert space, corresponding to 
an entangled state, is an open subset of an algebraic variety built by classical geometric constructions (tangent 
lines, secant lines) from the set of separable states. In this setting we describe well-known classifications 
of multipartite entanglement such as 2 x 2 x (n + 1), for n > 1, quantum systems and a new description 
with the 2x3x3 quantum system. Our results complete the approach of A. Miyake and makes stronger 
connections with recent work of algebraic geometers. Moreover for the quantum systems detailed in this 
paper we propose an algorithm, based on the classical theory of invariants, to decide to which subvariety of 
the Hilbert space a given state belongs. 

1 Introduction 

Let H = (g) C"^ ■ ■ ■ (E> C"*" be the Hilbert space of states of k particles. Denote by a basis of C"' with 
< ji < Hi — 1. Any state \'^) E H can be written as 

l<i<k Q<ji<ni-1 

where (g) • • • (g) is the standard basis of H. That basis will be denoted latter on by \ji . . .jk)- The 
Hilbert space % is an irreducibe representation (for its natural action [13]) of the semi-simple Lie group G = 
5i(ni,C) x ••• X SL{nk,'C). In the framework of Quantum Information Theory (QIT), G is the group of 
reversible stochastic local quantum operations assisted by classical communication (SLOCC, see |38|). and two 
states will be considered as SLOCC equivalent if they are interconvertible by the action of G, 

l^-) -^SLOCC I*) ^ I*) =g|$),|$) =g-i|*), with g = (51, . . . ,5^) e SL{ni,C) x ••• x SL{nk,C) 

Nonzero scalar multiplication has no incidence on a state l^*), therefore we can consider states as points in 
the projective space P(H) and SLOCC equivalent states will correspond to points in the same G-orbit. The 
representation "H of G has a unique highest weight vector which can be chosen to be w = |0 . . . 0) (it corresponds 
to a choice of orientation for the weight lattice see [IS|)- The orbit G.v C H is the unique closed orbit for the 
action of G on H and it defines a smooth algebraic variety after projectivization X = P(G.w) C P(?^). This 
variety X is known as the Segre embedding of the product of the projective spaces P"»~i, and is the image of 
the map (see [15]): 

(j): P(C"i) X P(C"^) x •■• X P(C"^) ^ P(C"i «)C"^ ® •••(8)C"^) 

(K], [^'2], ■ • ■ , [Wfc]) ^ [vi®V2«i---«iVk] 
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where Vi is a vector of C"' and [vi] the correspondhig point in P"' ^ = P(C"*). The variety X — ¥{G.v) = 
<?!)(P(C"i) X P(C"2) X • • • X P(C"^)) will be denoted by 

X ^ P"i"i X ■ • • X P""-! C ¥(n) 

From the QIT point of view [HIIS], the variety X is the set of separable states in F{H). 

In this paper we cover some examples (2 x 2 x (n + 1), with n > I, and 2x3x3 quantmn systems) of 
classifications of multipartite entanglement. We describe the entangled states by auxiliary varieties, i.e. varieties 
obtained from X by geometric constructions and we propose an algorithm to distinguish between the different 
states. The constructions of auxiliary varieties by secant and tangent lines are explained in Section [5] and some 
classical results of algebraic geometry on the dimension of those varieties are recalled. In Section [3] we give the 
geometric descriptions of the orbit closures for the quantum systems 2x2x(n+l),n>l, (Theorem [5]) and 
2x3x3 (Theorem [3]). The description of the (projectivized) Hilbert space by different classes of entanglement 
corresponds to a stratification of the ambient space by algebraic varieties with natural geometric inclusions 
among those varieties (the "onion like" structure of [W). Those inclusions are detailed in Figures [1] [51 Eland 
SI The proofs of Theorem |21 and Theorem |31 are based on dimension counts to identify the algebraic (auxiliary) 
varieties with the corresponding orbit closures. Those orbits are known from their representatives [31]. A 
geometric description of 2 x 2 x (n + 1), with n > I, quantum systems was already established in [391 using 
the concepts of projective duality, hyperdeterminants and singular locus of hyperdeterminants. In Section [31 we 
recover the results of [3S] directly from Theorem [51 We prove that Theorem [21 and the results of [22] are dual to 
each other and thanks to Theorem [51 we obtain more details in the geometric description of the singular locus 
of the dual variety. In Section [5l we introduce classical invariant theory in the context of hypermatrices. We 
obtain an algorithm to decide to which strata (variety) of the ambient space belongs a given state In other 
words our method allows us to identify the orbit of a given state j^*). This section is mainly based on classical 
invariant theory techniques but we relate part of the information obtained by those techniques with geometric 
descriptions. Finally we mention recent works of algebraic geometers which we believe should help to provide 
deeper understanding of multipartite entanglement. 

Notations 

We work throughout with algebraic varieties over the field C of complex numbers. In particular we denote by 

V a complex vector space of dimension iV + 1 and X" C P{V) — is a complex projective nondegenerate 
variety (i.e. not contained in a hyperplane) of dimension n. Given x a smooth point of X, we denote by T^X 
the intrinsic tangent space, T^X the embedded tangent space, of X at x (see [5T]). The notation X C V (resp. 
T^X) will denote the cone over X (resp. over T^X) and [v] G ¥{V) will denote the projectivization of a vector 

V G V. The dimension of the variety, dim(X), is the dimension of the tangent space at a smooth point. We say 
X £ X is a general point of X in the sense of the Zariski topology. The locus of smooth points of X is denoted 
by Xsrnooth and the locus of singular points by Xsing. 

2 Join varieties 

Let X and y be two points of P^, the secant line P^^ is the unique line in P^ containing x and y. We define the 
join of two varieties X and Y to be the (Zariski) closure of the union of the secant lines with x G X and y gY: 

j{x,Y)= y pi^ 

Suppose Y C X and let yy^ denote the union of P],'s where Pj is the limit of P^^^ with x € X , y € Y and 
x,y ^ yo G Y . The union of the ^ is defined as the variety of relative tangent stars of X with respect to 
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Y (see [S3]): 

T{Y,X)^ y Tly^^j 

yeY 

The following result due to F. Zak ([53J Chapter I Theorem 1.4) will be useful to analyse the stratification of 
the ambient space by auxiliary varieties: 

Theorem 1. Any arbitrary irreducible subvariety Y"' C X"', n > satisfies one of the following two conditions: 

a) dim{J{X, Y)) = n + m + 1 and dim{T{X, Y)) — n -\- m; 

b) JiX,Y)^T{X,Y). 

Remark 2.1. The expected dimension of J{X, Y) is rt + to + 1: there are n degree of freedom for the choice 
of y e y, TO degree of freedom for the choice oi x £ X and 1 degree of freedom for the choice of the point on 
Pl-y Therefore the previous theorem says that if J{X,Y) has the expected dimension then T{X,Y) has also 
the expected dimension and is distinct from J{X, Y). We will see that dimension calculations will be used later 
to prove the existence of certain statras. 

The dimension of J{X, Y) will be calculated with the following famous lemma (see [H] Chapter III). 

Lemma 1. [Teracini's Lemma] If z £ J{X,Y)smooth with z = x + y such that x £ Xsmooth, y e Ysmooth, 
then 

f,J{X,Y)=f,X + fyY 

Intersting particular cases arise when Y = X. The s-secant variety of a projective variety X C is the 
variety (Ts{X) defined to be the closure of the union of the linear span of s-tuples points of X 

a; 1 , . . . , a: s tE X 

where P^"^^;^ is a projective space of dimension s — 1 passing through xi,. . . ,Xs- In other words as{X) = 
J{X, as-i{X)) with (Ti{X) = X. The variety a2{X) is often called the secant variety. For the variety of relative 
tangent stars we obtain for Y = X the usual tangential variety: 

t{X)^T{X,X)^ U f,X 

Example 2.1. Let H = C'"+i ® C"+\ and X = P™ x P" c P(7^) be the Segre product of two projective 
spaces, i.e. the unique closed orbit for the action of SLm+i x SLn+i on P(H). The variety X corresponds to the 
projectivization of rank one matrices in the projectivization of the space of matrices of size (to + 1) x (n + 1). 
Therefore the s-secant variety <Js{X) is the projectivization of the set of rank less than s matrices (sum of 
s matrices of rank one). As noticed in [TB] the stratification of P(?^) by secant varieties of P™ x P" is the 
stratification by local rank for bipartite quantum systems : 

P™ X P" C (72 (P™ X P" ) C • • • C a™„(„,„)_i(P™ X P") C P(H) 

Remark 2.2. As noticed in [2], a projective line Pl.y in the Hilbert space P(?^) represents all possible super- 
positions of the states x, y £ T-L. 

Definition 2.1. Let X C P(y) be an irreducible variety of dimension n. The s-secant variety as{X) is said 
to be nondefective if either dim{as{X)) = sn + s — 1 or as{X) — P{V). 

A direct consequence of Theorem [1] is the following proposition : 
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Proposition 2.1. Let X" C f{V) he a nondegenerate variety. Let us assume that the k-th secant variety is 
nondefective and does not fill the ambient space. Then we have the following filtration with the given dimensions: 

^ (1t{X) Ca^iX) CT{X,a2iX)) CaaiX) (lT{X,a3{X)) C--- C ak{X) c¥{V) 

dim=7i ^2n =2n+l =3n+l =3n+2 =4n+2 =kn+k-l 

Proof. If ak{X) is nondefective, i.e. is of dimension kn + k — 1 then one knows from Theorem [T] that 
T{X,ak-i{X)) is of dimension kn + k — 2. Moreover Theorem [T] ensures us in this case dim{ak{X)) = 
dim((Tfc_i (X)) + dim{X) + 1, thus dim{ak-i{X)) = [k — l)n + (fc — 1) — 1, i.e. ak-i{X) is nondefective 
and we apply the same argument inductively □. 

Remark 2.3. Proposition 12. II gives a priori a filtration by secant and tangential varieties and explain part of 
the "onion like" structure described in [2>T\. As X is a Segre product of projective spaces, Definition 12.21 shows 
that there are others intermediate auxiliary varieties which will appear in the filtration of the ambient space. 

We will need the following definition to identify subvarieties of (72 (X) when X is a Segre product of irreducible 
varieties. In Section l4l Definition 12.21 will also be used to describe the singular locus of the dual variety. 

Definition 2.2. Let Yi C P"', with 1 < i < m be m nondegenerate varieties and let us consider X = Yi XY2X 
• • ■ X Y„i C p("i+i)("2+i)---(n,Ti+i)-i ^/jg corresponding Segre product. For J = {ji, . . . ,jk} C {!,..., m}, a J -pair 
of points of X will be a pair {x,y) G X x X .such that x — [vi(E)V2(^- ■ 'SJi'ji+i (Ei- ■ •(!S)Vj^(E)- ■ •(E)Vjf, ® • • -(SJWm] 

and y — \w\ ®W2® - ■ - ®Vj^ ® ''^ii+i Vj^ ® ■ ■ ■ ® Wr,i\, i.e. the tensors x and y have the same 

components for the indices in J . 

The J -suhsecant variety of a2{X) denoted by 02^X1 x ■ ■ ■ x Y_j^ x ■ ■ ■ x F^-^ x ■ ■ ■ x Y^) x Yj^ ^ ^2 ^ ' ' ' ^ ^ & 
is the closure of the union of line ^^.y with (x, y) a J -pair of point: 



f 2(^1 X • • • X Y_j^ X ■ ■ ■ X Y_j^ X ■ ■ ■ X Ym) X Yj-^ X Yj^ X ■ ■ ■ X Yj^ — ^(x,y)eXxX,(x,v)J-vair of points^xy 

Remark 2.4. The underlined varieties in the notation of the J-subsecant varieties correspond to the common 
components for the points which define a J-pair. Roughly speaking those components are the "common factor" of 
X and y in the decomposition oi z = x-\-y £ (T2 (5^1 x • • • x x • • • x Y-j^. x • • • x Ym )xYj^x Yj^ x - ■ - xYj^,. For instance 
when we consider the {l}-subsecant (respectively the {m}-subsecant) variety we can indeed factorize the first 
(respectively the last) component and we have the equality cr2(lli x F2 x ■ • • x Ym) xYi —YiX 0-2(^2 x • ■ • x Ym). 

Remark 2.5. For J = 0, the J-subsecant variety is a2{X). 

3 Stratification of the multipartite entangled states by the tangential, 
secant and join varieties 

We now state with Theorem [5] and Theorem [3] our geometric descriptions of entangled states for quantum 
systems of type 2 x 2 x (n + 1), n > 1, and 2 x 3 x 3 by join, secant and tangential varieties. The orbit closures 
will be denoted where the subscript is either a representative or a roman number to identify the orbit. 

Theorem 2. For quantum systems of type 2 x 2 x (n + 1), there are 6 (n = 1), 8 (n = 2) and 9 (n > 3) SLOCC 
entangled classes. Each entangled state corresponds to an open subspace (smooth points) of an algebraic variety 
build up from X = x x P", the unique closed orbit for the action of G = SL2 x SL2 x SL^+i on 
P(7^) — P(C^ ®C^~^^) by join and tangential varieties. The identifications of those algebraic varieties are 
given in Tables\^ [H and[^ and the partial order among them is represented in Figures\^ [H [7} 
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Orbit closure 


Normal form (representative) 


Variety 


Dimension 


Ovi 


|000) + |111) 


p7 


7 


Ov 


|100) + |010) + |001) 


T(P1 X Pl X Pl) 


6 


Oiv 


|001) + |111) 


(72 (P^ X Pl) X Pl 


4 


Oiii 


|100) + |111) 


Pl X Cr2(P^ X Pl) 


4 


On 


010) + |111) 


0-2(1°^ X F X pi) X pi 


4 


Oi 


|000) 


pi X pi X pi 


3 



Table 1: Identification of orbit closures and varieties for the 2x2x2 quatum system 



Orbit closure 



Normal form (representative) 



Variety 



Dimension 



Oyill 
Oyil 
Oyi 
Ov 

Oiv 
Oiii 
Oil 
Oi 



|000) + |011) + |101) + |112) 
|000) + |01f) + |102) 

|000) + |111) 
|000) + |0I1) + |101) 
|000) + |011) 
jooo) + jlOl) 
jooo) + jllO) 
|000) 



pn 

J{X,Oiv) 
<y2{x) 

t{X) 

pi X (T2(P^ X P2) ~ Pl X P5 

(T2(P^ X F X P2) X Pl 
(72 (P^ X Pl) X P2 ~ P3 X P2 
X = Pl X X P2 



11 

10 
9 

8 
6 
6 
5 
4 



Table 2: Identification of orbit closures and varieties for the 2x2x3 quatum system 



Orbit closure 



Normal form (representative) 



Variety 



p4n+3 

J{X,Oiv) 

02{X) 

t{X) 

»1 X cr2(Pl X P") 2± P 

cr2(Pl X Pi X P") X Pl 
(72(Pl X Pl) X P" ~ P3 X ] 



Dimension 



Oviii 

Oyil 
Oyi 
Ov 

Oiv 
Oiii 
Oil 
Oi 



|000) + |011) + |102) + |113) 
|000) + |011) + |101) + |112) 
|000) + |011) + |102) 

|000) + |111) 
|000) + |011) + |101) 
|000) + |011) 
jooo) + jlOl) 
jooo) + jllO) 
|000) 



j2n+l 



4n ■ 

3n 

3n ■ 

2n 

2n 

2n 

2n 
nH 
n H 



-3 
-5 
-4 
-5 
^4 
^2 
-2 

3 

2 



Table 3: Identification of orbit closures and varieties for 2 x 2 x (n + 1) quatum systems (n > 3) 



Remark 3.1. After [T^] the classification of entangled states of the 2x2x2 quantum system received a lot 
of attention because it showed for the first time that three qubits could be entangled in different inequivalent 
states. The existence of the so-called GHZ-state and W-state for a 3-qubit system proved in [12J is equivalent 
in our description to the existence of the (nondefective) secant and tangential varieties of P^ x P^ x P^ which 
is a classical result of algebraic geometry. 

Remark 3.2. In [3] different geometric characterizations of the 3-qubit system are given. The authors consider 
the intersection of the variety of pure state X with the three plane of symmetric tensors (symmetric states). 
That intersection reduces to the twisted cubic C C Psy,„ and allows them to characterize the GHZ-states and 
W-state in terms of that curve. 

Remark 3.3. The geometric classification of entangled states for 2 x 2 x (n + 1), n > 2, systems was proposed 
in [37. .39j but with a different geometric perspective (i.e. by dual varieties). Not all stratas are geometrically 
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p7 



r(pi X Pi X pi) 




Figure 1: Stratification of tlie ambient space for the 2x2x2 quantum system 



pii 
J{X,Oiv) 




X = pi X pi X p2 



Figure 2: Stratification of the ambient space for the 2x2x3 quantum system 

described in the paper. We wih make the connection between our identifications and the classifications of [571155] 
in Section m 

The next theorem provides a geometrical description for 2x3x3 quantum system. 

Theorem 3. For a quantum system of type 2x3x3 there are 17 different SLOCC entangled classes. Each 
entangled state corresponds to an open subspace of an algebraic variety constructed from X = ]P^ xV'^ x P^, the 
unique closed orbit for the action of G = SL2 x SL3 x SL3 on P(H) = P(C^ Cg) (8) C^), by join and tangential 
varieties. The identifications of those algebraic varieties are given in Table and the partial order among them 
is represented in Figure^ 
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p4ri+3 



a3(X) 

JiX,Oiv) 




X = pi X pi X 



Figure 3: Stratification of the ambient space for 2 x 2 x (n + 1), n > 3, quantum system 

Remark 3.4. The geometric inclusions between the varieties of Tables [TJ [21 [3] and U are represented in Figures 
[U [51 [3]and[H Those inclusions can be deduced from the relations between the representatives of each orbit, 
e.g. in Table [H it is clear that Ojx is a subvariety of Oxiii- But it could also be deduced from geometric 
considerations, e.g. the tangential varieties are included in the secant varieties (i.e. Ov C Ovi)- The geometry 
of the auxiliary varieties reflects how inequivalent entangled states are partially ordered under local actions. 
This was emphasized in |37[ I39j and it becomes very natural in our description by auxiliary varieties. 

Remark 3.5. As we point it out in Section[Sla lot of classifications of entanglement in quantum systems can be 
deduced from older results from representation theory and invariant theory. For instance for the 2x3x3 system, 
the classification was first established in QIT context in [T^ but can be deduced from [IT] . It is interesting to 
notice that the classifications of tripartite entangled states of |10ll571[55l[55] can be obtained from a classification 
theorem on trilinear forms of type 2 x n x n proved by Camille Jordan in 1907 [22) . 

Proofs of Theorem [2] and Theorem [3] 

We now prove the identifications of the orbit closures with algebraic varieties constructed from the unique closed 
orbit under the action of G = SLp x SLg x SLr on P(Cp «) C? C) for 

{2x2x2 Table[Tl 

2x2x3 Tabled 
2 X 2 X (n + > 3 Table[31 

2x3x3 Table[l 
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Orbit closure 


Normal form (representative) 


Variety 


Dimension 


Oxvii 


|000) + |011) + |100) + |122) 


pl7 


17 


Oxvi 


000) + Oil) + 101) + 122) 


J{X,t{X)) 


16 


Oxv 


|000) + |011) + |022) + |101) + |112) 


T{X,t{X)) 

J(X,P1 X (72 (P^ X P2)) 


15 


Oxiv 


000) + Oil) + 122) 


14 


Oxiii 


|000) + |011) + |022) + |101) 


r(X,pl X cr2(P^ X P2)) 


13 


Oxii 


|000) + |011) + |101) + |112) 


CT2(0-2(P^ X F X P2) X P2) 
J(p5 X p2^ cr2(P^ X F X X P2) 


13 


Oxi 


000) + Oil) + 121) + 102) 


13 


Ox 


|000) + |011) + |102) 


J(X,cr2(P^ XF XP2) XP2)) 
Pl X cr3(p2 X P2) ~ Pl X P^ 


12 


Oix 


000) + Oil) + 022) 


9 


Oviii 


|000) + |011) + |110) + |121) 


(p5 X P2) 


13 


Ovii 


|000) + |011) + |120) 


J(X,P^ X P2) 


12 


Ovi 


|000) + |111) 


CJ2{X) 


11 


Ov 


|000) + |011) + |101) 


r{X) 

Pl X 0-2(P^ X P2) 


10 


Oiv 


|000) + |011) 


8 


OlH 


000) + 101) 


(T2(P^ X F X P2) X P2 


7 


On 


000) + 110) 


cr2(P^ X P2) X P2 - P5 X P2 


7 


Oi 


|000) 


X = Pl X P2 X P2 


5 



Table 4: Identification of orbit closures and varieties for 2 x 3 x 3 quatum system 



The normal forms (representatives) of the orbits under the action of G are given in [41] and there is a finite 
number of them. We express those normal forms in the tables with the braket notation by the convention of 
the introduction \ijk) = ei ® ej (E> where (et)o<t<s-i stands for a basis of C*. The first step in our proof is 
to calculate the dimension of each orbit closure. It can be done directly using the following two facts: 

• Let Ox — P{G.x) C F{V) be the closure of a G-orbit with representative x, for G a semi-simple Lie group 
with Lie algebra q. Then T^O = q.x. The notation q.x represents the action of g on the representation V 
(|13|). Moreover ii x = u + v one has T^O = g.x = q.u + q.v. 

• Consider the Segre product X = P(G'.(ei ® ■ • • (g) e™)) = Fi x • • • x y„ C P{Vi ® • • • (g) Vm) with G a 
semi-simple Lie group, then we have Te^(g)...ge„^ = Jei^L 62 ® ■ • ■ €5 em + ei eg) Te^Y2 g) • ■ ■ g) + 
h ei (g 62 (g • • • g) Te^Ym- In particular for X = PP^^ x P^^^ x P''^^ and e (g / (g g £ X, we have 

= CP g) / g) g -t- e g) C« g) 5 -I- e g) / g) C. 

Those remarks allow us to calculate the dimension for any orbit from the normal form. For instance for the 
orbit closure Ovi we get: 

71ooo)+|iii>G.(|000) + |111)) = 7]ooo)G.|000) + f|ni)G.|lll) 

= g) ep g) ep + ep g) g) ep -f ep gi ep g) C'' + g) ei g) ei + ei g) g) ei ei (g ei g) C''^ (1) 

dim=2(p+(}+r)-4 

After projectivization dim(G.(|000) + |111))) = 2{p + q + t) — 5. That orbit is clearly a subvariety of (T2(PP^^ x 
pg-i X p'-i) and dim(g2(^)) < 2 (p - l + q- l + r-l) + l. The variety cr2(X) is irreducible because X 
is. Thus we have G.(|000) + |111)) C o'2(P'^^^ x P'*'^^ x P*"^^) with equahty of dimensions which proves the 
equality G.(|000) -I- |111)) — o'2(PP^^ x P''^^ x P''^^). The secant variety is of maximal dimension then one 
deduces from Proposition 12.11 that the tangential variety t{PP^^ x P''^^ x P*"^^) is of dimension 2{p + q + r) —6 
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pl7 

J{X,t{X)) 




X = Pl X P2 X P2 

Figure 4: Stratification of the ambient space for the 2x3x3 quantum system 



and we identify it with the orbit Oy which is the only one of dimension 2(p + q -\- r) — For the orbits On 
and Oiv one notices that |000) + |011) = eo (eo (8> eo + ei ei) and |000) + |110) = (eo eo + ei (8> ei) eo 
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which allows us to identify those orbit closures with cr2(P^"^ ® P"^"^) x P'^-^ and PP-^ x cr2(P'^"^ x P'-^) after 
a dimension count (again the orbit is clearly a subvariety of the {l}-subsecant or the {3}-subsecant and has the 
same dimension). The orbit closure P(G.(|010) + |111))) is the orbit of the line Pe(,(8)ei(g)eo,ei(giei®ei defined by 
the {2}-pair of points (cq (8) ei (g) eo, ei ei ® ei). The corresponding orbit closure is the {2}-subsecant variety 
CT2(P^ X Pi X Pi) X Pi. This completes the identifications of Table [T] _ 
For Tables [2] and [3] we observe, after a dimension calculation similar to ([T]), that the dimension of Ovii is 
equal to 3ri + 4 [n — 2 for Table [5] and n > 3 for Table [3]) . Moreover it is clear from the representatives that 
J{X,Oiv) 3 Ovii- That last assertion comes from 

XVII + |011) + |102) e J{X,Oiv) 

We also know that dim( J(X, Oiv)) < dim{X) + dim(0/y) + 1 = 3n + 5. Thus we have Ovii C J{X, Oiv) and 
3n + 4 = dim(Oy//) < dim( J(X, Oiv)) < 3n + 5. We now prove with Teracini's Lemma that J(X, Oiv) is in 
fact of dimension 3n + 4. Let x ^ e ® f ® g -\- h ® {m ® n ~\- p ® q) be a smooth point of J{X, Oiv)- Teracini's 
Lemma says that the tangent space of the join is given by TxJ{X, Oiv) — Te®f®gX + T\i^{^jn®n-\-p®q)0 iv with 
fe^f^gX = C2(K)/(8)5r+e(X)C2(K)5r+e(K)/(X)C"+i &ndfh^(^rn(g,n+p(g,q)Oiv = C2®(m®n+p(8)g) + /i®C2"+2. When 
we look at the intersection of the tangent spaces we have h ® f ® g ^ Te^f,g,gX n Tj^^(^jyi^n+p®q)0 iv ^ thus the 
intersection does not reduce to {0} and therefore the join is not of maximal dimension, i.e. dim( J(X, Oiv)) l£ 
3n + 4. We conclude that dim(J(X, Oiv)) = 3n + 4 and corresponds to Ovii- The orbit Oviii is of dimension 
11 for Table [5] and corresponds to the ambient space P^^. The orbit Oviii is of dimension 3n + 5 for Table [3] 
but it is known, [7J, that 0-3 (P^^^ x P"?^^ x P''"^) is of dimension 3p + 3q + 3r — 7 and therefore we can state 
that Oviii = asiP^ x P^ x P"). 

The last identifications concern Table HJ with Teracini's Lemma one shows that J{X,t{X)) has the expected 
dimension, i.e. 16. Then the orbit closure Oxvi, which is the unique orbit of dimension 16, corresponds to 
J{X,t{X)) (indeed we have J{X,t{X)) = T{X,a2{X)) because here a3{X) fills the ambient space). One 
deduces there exists an orbit closure of dimension 15 which corresponds to T{X,t{X)). This has to be the 
orbit closure Oxv (the only orbit of dimension 15). It is clear from the normal form that Oxiv is included in 
J(X,pi X cr2(P^ x p2)). But the expected dimension of J{X,¥^ x (T2(P^ x P^)) is 14 which is the dimension of 
Oxiv- Thus we conclude to the equality between Ox/v and J{X,¥^ xa2{P'^ xP'^)). As J(X,pixcr2(P^xp2)) has 
the expected dimension one knows there exists a variety of dimension 13 corresponding to T{X, P^ x (T2 (P^ x P^)) . 
But according to [H] there is only one orbit of dimension 13 in Oxiv and that is Oxiii- Thus the orbit closure 
Oxiii corresponds to T(X,pi x cr2(P^ x P^)). The varieties cr2(Pi x P^) x P^ and cr2(F^ x F x P^) x P^ are 
isomorphic to P^ x P^ and 172 (P^ x P^) is of dimension 13 (that's the projectivization of the set of rank at least 
2 matrices in the projectivization of the space of 6 x 3 matrices). But there are three others orbits of dimension 
13 which are orbits Oviii-, Oxi and Oxii- From the normal forms we can affirm that orbit Oxii is contained 
in CT2(ct2(P^ X Pi X P2) X P2) and the orbit Oviii is contained in cr2(cr2(P^ x P^) x P^) = cr2(P^ x P^) because 

xx// = |000) + |011) + |101) + |112) = |000) + |101) + |011) + |112) 
XVIII = |000) + |011) + |110) + |121) = |000) + |011) + |110) + |121) 

e<T2(PlxP2)xP2 ecr2(PlxxP2)xP2 

That leads to the equalities Oxii = cr2(CT2(P^ x F x P^) x P^) and Oy/// = cr2(P'"' x P^). The last orbit closure 
of dimension 13, namely Oxi, is a subvariety of J(P^ x p2,cr2(P^ x F x P^) x P^): 

XXI = |000) + |102) +|011) + |121) 

e(T(pixPixp2)xp2 ep5xp2 
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By Teracini's Lemma one obtains that J(P^ x P^, (T2(P^ x x P^) x P^) is of dimension less than 13. Indeed let 
M=(a(g)6 + c(8)fi)i8)eeP^xP2, then fuP'^ x = C*^ O e + (a (g) 6 + c (g) d) (g) and v = f ® g®h + k® g®l & 
cr2(IP^ x F X P2) x P2, then f„CT2(P^ x f^yi x = ® g ® C'^ + {f ® g® h + k ® g® I) ® . Thus 
(f„p5 xp2)n(T„cr2(F^ xIP^xp2) xp2) D C2(g)g(g)e, i.e. the dimension of the join variety is at most 15-2 = 13 
and therefore we have Oxi = J(F^ x P^, CT2(Pi x F x P^) x P^). 

The same argument holds for orbits Ovii and Ox- One first shows that Oyii C J{X,f"^ x P-^) and Ox C 
J(X,cr2(Fi x F X P2) X P2): 




GCT2(PlxP2)xP2 ex 



XX = |000) + |102) +|011) 

e(T(pixPixp2)xp2 gx 

Then Teracini's Lemma allows us to prove that the join varieties are of dimension less than 12 and we conclude 
to the equality. 

The last orbit to identify is Ojx which is of dimension 9. From the normal form this orbit is clearly included 
in P^ X cr3(P^ X P^) = P^ X P^ and the identification follows because of the equality of dimensions. □ 

4 Back to Miyake's geometric description by dual varieties 

Since G acts with a finite number of orbits on P('H), the orbit structure of P(H*) is identical to the orbit structure 
of P('H). In this section we identify the orbit closures in P(H*) with duals of varieties of the stratification given 
by Theorem We then recover Miyake's geometric descriptions of entangled states by dual varieties for the 
2x2x2 and 2 x 2 x (n + 1) quantum systems. 

The dual variety and its singular locus 

Let X C P(y) be a projective variety and let T^X denote the embedded tangent space of X at a;, a smooth 
point of X . Define the dual variety X* by 

X* = {H e P(y*)| 3xe Xsmooth such that f^X C H} C ¥{¥*) 

The biduality theorem {X*)* — X, true in charestic zero, implies that the original variety can be reconstructed 
from its dual variety. The dual varieties have been studied intensively by algebraic geometers (|46j). In the 
case of the variety X = P" x P" c p("+i) ~^ (projectivization of rank one matrices) it is well known that the 
dual variety can be identified with the variety of matrices of rank at most ni. Thus, up to multiplication by a 
nonzero scalar the equation defining X* is the determinant. This leads to a higher dimensional generalization of 
the determinant, called hyperdeterminant, which was first introduced by Cay ley and rediscovered by Gelfand, 
Kapranov and Zelevinsky [T3]. The hyperdeterminant in the sense of [13] is the defining equation of the dual 
of X = P"i x • • • X P"*" when X* is a hypersurface. 

In |[37j A. Miyake uses this notion of hyperdeterminant to classify multipartite entangled states for 2 x 2 x n 
quantum systems. In the bipartite case, the dual of the set of separable states, P" x P" c p("+i)'-i, is 
isomorphic to the the projectivization of the set of matrices of rank less than n. The generic entangled state, 
in the bipartite case, corresponds to matrices of maximal rank and therefore corresponds to p("+i) 
Following the analogy with the 2-dimensional case Miyake proposes to use the stratification of P{V*) by X* 
and its subvarieties to distinguish the states of entanglement. The variety of separable states X being SLOCC 
invariant so is X* and its singular locus. The dual variety and its singularities induce a filtration of the (dual) 
ambient space: 
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X*„^ CX*C P{V*) 

In order to explain what the singular locus of the dual variety is, we need to look at the tangent hyperplanes of 
the variety X. When X* is a hypersurface, a smooth hyperplane H £ X* is a hyperplane tangent to X with 
a unique singular point which is a nondegenerate quadric. In other words the restriction to X of the linear 
form which defines _ff is a quadric of full rank (a hypersurface with a Ai singular point). Thus when X* is 
a hypersuface there are two ways for a hyperplane to not be a smooth point, either by having more than one 
point of tangency or by defining a degenerate quadric |20[ (55] . 

Definition 4.1. Let X C ¥{V) and X* be its dual variety, which we assume to be a hypersurface. We define 
-^node' node component of X* to be the set of hyperplanes having more than one point of tangency : 

Kode = {HeX*, 3ix, y)eXxX,x^y,fxXc H, fyX C H} 

We define the cusp component X*^^^, to be the set of hyperplanes defining a singular hyperplane section with 
degenerate quadratic part : 

Kups = {H eX*,3xGX, f^X cH,{Xn H, x) ^ A^} 
Following [14] we can decompose X'^^j^^ into irreducible components : 

Definition 4.2. Let X = Y"^^ x . . . F"™ C p{"i+i)---("m+i)-i he the Segre product ofm nondegenerate varieties 
and J = {ji,...,jk} C {1,. . 

J) = {H e ViV*),3ix, y) aJ^ pair of point of X x X, f.,X C H, f.,X C H} 

We now prove a proposition which describes for Segre products the node components in terms of the dual of 
the subsecant varieties. 

Proposition 4.1. Let X = Yi x • • • x Ym C P(\^) be the Segre product of m nondegenerate varieties and let X* 
be its dual variety. The J -node component is the dual of the J -subsecant variety, i.e. for J — {ji, . . . C 
{1, . . . ,m}.- 

Kode{J) = (^2(^1 X • • • X y^, X • • . X y^.^ X • • . X y„o x Y,, X y,, X • . • X y,j* 

Proof. It is a consequence of Teracini's Lemma. Let us denote by Z the J-subsecant variety Z = a2{Yi x • • • x 
Hji X • • • X y^^ X • • • X Yjn) X Yj-^ X y,2 x ■ ■ • X Yj^ and z^x + y^Za general point of Z. By definition {x, y) is 
a J-pair of point. According to the Teracini's Lemma we have T^X + TyX — T^Z . Thus if iJ e Z* is tangent 
to Z at z it means H D T^X and H D TyX i.e. H £ X'^^^^{J) because [x, y) is a J-pair of point. On the other 
hand if iJ G X*^^^{J) then there exists a J-pair of point {x,y) such that H D T^X and H D TyX i.e. H is 
tangent to Z at the point z ^ x -\- y, i.e. H E Z* . □. 

Remark 4.1. This proposition shows in particular that X*^^^{(d) = a-2{X)*. This equality is used in [20] to 
study the dimension of the singular locus of the duals of Grassmannians. 

Stratification by the dual variety and its singular locus 

We now recover Miyake's classifications j37[ [55] by dual varieties by establishing the isomorphisms between the 
varieties of Theorem [51 and the varieties in the dual space. Morevover the components of the singular locus of 
the dual variety are described in terms of the duals of auxiliary varieties. 

Theorem 4. For 2 x 2 x (n -f- 1) (n > 1) quantum systems the duality between the orbit closures are given in 
Tables\M\Eand\l\ 
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Orbits 


Varieties 






r(pi X F 


1 X pi) 2± (pi X F 


1 X F 




Oiv ^ 0]y 


(CT2(P1 Xl 


ji) X pi) 


~ ((T2(P^ X pi) X 


pl)* 




Oui ^ 0]jj 


(Pl X (72(1 


Jl X pi)) 


~ (Pl X cr2(P^ X ] 


pi))* 


= ^:ode({i}) 


On - 0*jj 


(cr2(P^ X Pi X 1 


ji) X pi) 


~ (cr2(P^ X Pi X ] 


pl) X 





Table 5: Duality between orbit closures for the 2x2x2 quantum system 



Orbits 


Varieties 


OvH ^ O] 


J(pi X pl 


X P2,I 


pl X (7(P1 X P2)) 


- (pl X pl X p2)* 


Ovi ^ 0*u 


(72 (P^ X pl X p2) : 


^(^2( 


pl X pl) X F2)* - 


(P3XP2)*=X*„,,({3}) 




r( 


pl X P 


1 X P2) ~ T(P1 X 


pl X p2)* 


Oiv ^ 0]y 


pl X (72 (pl 


X F2) 


2± (pl X cr2(Pl X I 


^'))* = ^,:o.e({i}) 


Oui ^ 0\u 


(72 (pl X Pi X P2) 


X pl 


~ ((72(Pl X Pi X I 


^2) X pl)* = x,:„,,({2}) 



Table 6: Duality between the orbit closures for the 2x2x3 quatum system 



Orbits 


Varieties 


Oviii — Oil 


(73 (Pl 


X pl X P") ~ ((72(Pl X pl) X P")* ~ (P3 X P")* 


Ovn ^ 0*i 


J(Fi 


X pl X P",Fl X (72(Pl X P")) ~ (Pl X pl X P")* 


Ovi ^ Ovi 




C72(Pl X pl X P") ~ (72(Pl X pl X P")* 


Ov^Ov 




T(Pl X pl X P") ~ T(P1 X pl X P")* 


Oiv - 0*iv 


(Pl X (72 (Pl X 


P") ~ pl X P2"+1 ~ (Pl X p2"+l)* ~ pl X ((72(Pl X P"))* 


Oui ^ 0*iii 


(72(1 


='1 X Pi X P") X pl ~ ((72(Pl X Pi X P") X pl)* 



Table 7: Duality between the orbit closures for the 2 x 2 x (ti + 1) quantum systems, n > 3 



Proof. Most of the identifications follow from calculation of the dimension of the dual of each variety of Theorem 

121 

1. Table [SJ it is well known ([5S]) the dual of pi x pi x pi is a hypersuface and therefore corresponds to the 
closure of the unique orbit of dimension 6, i.e. (pi x pi x pi)* ~ T(pi x pi x pi). 

The {j}-subsecant varieties are isomorphic to P^ x pi. But P'^ x pi is the projectivization of rank one 
4x2 matrices. The set of rank one 4x2 matrices is equal to the set of degenerate 4x2 matrices (the 
generic 4x2 matrices are of rank 2 and the degenerate ones of rank strictly less than 2). It means P"^ x pi 
is self dual, i.e. P'^ x pi ~ (P'^ x pi)* and therefore so are the {j}-subsecant varieties. 

2. Table [B] here again the dual of pi x pi x P^ is a hypersurface and thus according to Theorem [5] one gets 
(pi X pl X p2)* ~ J(pi X pl X p2,pi X (72(Pi X p2)) (the uuiquc hypersurface). The orbit closure On 
is isomorphic to P^ x P^. But the projectivization of the set of rank one 4x3 matrices is dual to the 
(projectivization of the ) set of rank at most 2 matrices which is of dimension 9 (after projectivization). 
Thus we identify On and the unique orbit closure of dimension 9, which is Oyi- The orbit closures Om 
and Oiv are self-dual because isomorphic to pi x P^. The remaining variety Oy has to be self-dual. 

3. Table [71 in this table the dual of pi x pi x P" {n > 3) is no longer a hypersurface. A dimension count, 
using Katz's dimension formula ([Ml), shows that (pi x pi x P")* is of dimension 3n-t-4 and therefore is 
isomorphic to J(pi x pi x P",pi x (72(Pi x P")) according to Theorem[21 The orbit On is isomorphic to 
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p3 ^ pn ^jjg projectivization of rank one 4 x (n + l) matrices. Its dual variety is the projectivization of the 
4 X (n + l) matrices of rank less than 3. But the projectivization of 4 x (n + l) matrices of rank less than 
3 is of dimension 3n + 5 and therefore Ojj ~ Oviii- The orbit closures Om and Ojv are isomorphic to 
pi ^ p2n+i therefore are self dual. To conclude one calculates the dimension of 0-2(1?^ x x P")* using 
techniques proposed in [50] where the dimension of a2{X)* has been studied. Those calculations lead to 
dim(cr2(P^ X P^ X P")*) = 2n + 5 and therefore the secant variety is self-dual. It forces the tangential 
variety to be also self-dual. □ 

Theorem fallows us to recover Miyake's classification of 2 x 2 x (n + 1) quantum systems for n > 1. For instance 
a direct consequence of Table [S] and Figure [T] is Figure [5] which is the geometric description developped in [?7] . 

(P7)* 



(pi X pi X pi)* 



(cr2(P^ X Pl) X Pl)* 

pi X pi X pi ~ r(pi X Pi X pi)* 

Figure 5: Stratification by the dual variety and its singular locus 

Moreover with Proposition 14.11 the node components are identified. Those node components turn out to be 
self-dual varieties and are geometrically described. The same comment is true for the 2x2x3 system. For 
2 X 2 X (n + 1) with n > 4 the interpretation of the node components is less obvious as the dual variety of 
pi X pi X P" is not a hypersurface. 

5 Algorithms 

In this section, we describe an algorithmic method to identify the orbit of a given state l^*). Consider G a 
semi-simple Lie group and V a representation of G. We know from Kac's classification [23] which pairs (G, V) 
have a finite number of orbits. A method to classify those orbits is proposed by Vinberg ^50^ .Slj and leads to the 
determination of normal forms. Let Vi = €.'"■' , with 1 < i < k, he k complex vector spaces. The classification 
of the orbits is well known when G is a product of linear groups, G — GL{Vi) x • • • x GL{Vk), acting on the 
tensor space TL :— Vi ^ ■ ■ ■ ® Vk with finitely many orbits (or equivalently, SL{Vi) x • • • x SL{Vk) acting on 
the projective space P{TL)). This is the case when fc = 2 (it reduces to the classification of matrices according 
to their ranks, see example 12. ip and fc = 3 with ni x n2 x na = 2 x 2 x n and ni x n2 x na = 2 x 3 x n. The 
normal forms are given in [41J and in the previous sections we took advantage of the knowledge of those normal 
forms to describe the orbit closures. When G acts with an infinite number of orbits, some special cases where 
the normal forms depend on parameters are solved, see for instance the 3x3x3 case in [40 . 
That same problem of classification of orbits under the product of linear groups has been studied from the 
point of view of classical invariant theory up to about 1940 (see e.g. [T71 [151 [IHl 133 HZl |3S| ) and suscited quite 
recently a regain of interest, motivated by potential application to quantum computing (see e.g. [5l [6l I251I371I391 
138]). The classical invariant theory approach has the advantage to produce polynomials (invariants, covariants. 
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concomitants) which can be used to separate the orbits. For instance the case 2x2x2 goes back to 1881 with 
the work of Le Paige [21] where a complete hst of covariants is described. That hst of covariants allows us to 
distinguish the orbits. The 2x2x2 case was also treated into a more precise form by Elise Schwartz |43j and 
independently by Saddler [12] • The result is reproduced in Sokolov's book [U] where the state of the knowledge 
up to 1960 is summarized. The classification for the 3x3x3 case, following the same techniques, can be found 
in |47j . Very little is known about other cases in terms of invariants, covariants and concomittants. 
For our purpose a state 

l<i<k 0<ji<ni-l 

will be assimilated to the hypermatrix A = i^ji,...,jk)o<j <n -i vi=i k ^^"^ first covariant of the hypermatrix 
will be the multilinear form (also called the ground form): 

A{x,y, . . . , z) :— ^ ^ Aj^^^^^j^Xj^^yi^ ■ . ■ zj^ 

l<i<k 0<ji<ni-l 

In principle the classification of the orbits can be obtained from the knowledge of the invariants, covariants 
and other concomitants (in the sense of classical invariant theory) of hypermatrices. This is a non-trivial task, 
and the result would allow in particular to write down explicit equations for the orbits closures, i.e. for the 
corresponding algebraic varieties. 

We will mainly make use of the method of Schur functions, introduced in invariant theory by D.E. Littlewood 
[32]. Our strategy will be to find rigorously, or by a guess from numerical data, a generating function for the 
number of covariants of any given type (the number of general concomitants is more difficult to compute), and 
then, guided by the series, to try to construct explicitly first the covariant polynomials and hence some other 
concomitants by all possible methods. 

Once the concomitants are obtained, the classification of the orbit closures can be recovered by testing the 
nullity of the concomitants. The description of the algebra of concomitants is a very tedious process which 
needs (even for simple cases; for general case this is unrealizable) several hours of computations on a computer 
algebra system. Nevertheless, once the polynomials are obtained and written in an appropriate way, the test is 
very efficient since this amounts to evaluate the polynomials on representatives of the orbits. 

5.1 General method 

When it is possible, we first determine the Hilbert series. The Hilbert series is easier to obtain than the 
description of the algebra and it will allow us to guide the calculations. 

The sets of all invariants and covariants of hypermatrices of a given size are algebras Inv = S{H)^ and Gov = 
[S{H) dg) S{V{ © • ■ • © V"^*)]*^ which can be graded according to the degree d. The action of G on the space Gov 
provides an additional information: the weight of a covariant regarded as a relative invariant. More precisely, 
the weight of an invariant F G S'^{'H) is the vector £ — {£i, . . . , £k) such that for any g — {gi, . . . , gk) G G 

g.F^idetgiY' ...{detgkY''F, 

where g.F means the image of F under the natural representation of G on S'^{'H). Similarly a covariant of 
degree d — {dQ,di, . . . ,dk) is a relative invariant of S{TL) ^ S{V{ © • • • © V^*). So the algebra of the invariants 
and covariants can be graded according to both the degree d and the weight £. 

For simplicity we will consider only the space Inv(cio) (resp. Cov((i)) of the (resp. multi) homogeneous polyno- 
mials of fixed degree do (resp. multi-degree d) and we study the multivariate Hilbert series 

Hcov{t;u) := ^dimCov(d)i'^°Mf^..u;?\ 

We show (see appendix [2] for details) the Hilbert series can be written using a Cauchy function: 

n.I^1-n(T^)°" 

rn. ^ 
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where S — oi„im is a (potentially infinite) linear combination of certain elements m. For instance 



3 , 1 - tu-^ 



2 {l-tvY{l-tuv)i 
We have: 

Hcov{t;u) = CT„jn>' ...CTi,„ri>"i?i^(ui,wi)...i?i„(u„,w„)nt •• -Ai^ (■«„,«„)] (2) 

where CTyf{v) denotes the constant terms of the Laurent series f{v), fly is the Macmahon operator which 
sends the negative power of v to 0, A2{u,v) — u + ^, B2{u, u) = 1 — for binary variables and 

A2iu,v) = u + v + B3{u,v) = (l - J-) (l- J_) (l-- 
vu \ vw^ ) \ v^u J \ u 

for ternary variables. 

Note using equation (fT4)) to compute a closed form for the Hilbert series is not straightforward. The main strat- 
egy consists in decomposing the Laurent series into simple fractions, sending the fractions which contributes 
negative powers of each u to 0. This is a very tedious calculation which can be performed only for the simplest 
cases. 

Now, our method of construction consists in generating concomitants and testing if the dimensions of the graded 
spaces spanned by these coincide with the dimensions predicted by the Hilbert series. 

Once constructed, we use the covariants to identify the orbit of a given form. Indeed, if P is a concomitant 
then the assertion P = is invariant on an orbit. So the goal is to compute sufficiently many concomitants in 
order to distinguish the orbits. Since in our case the classification is already known, we do not need to describe 
completely the algebra. 



5.2 The case 2x2x2 



To illustrate the method, we first apply it to the (well known) simplest non-trivial case: fc = 3 and Vi — . 
The generating series of the algebra of covariants is known : 



t^u\u\u\ 



(1 - toiU2-U3)(l - f^ul){l ~ f^uDiX ~ t'^ul){\ - t^UxU2U-i){\ - t^) ' 

This suggests that the algebra is generated by a trilinear covariant of degree 1, three quadratic covariants of 
degree 2, a trilinear covariant of degree 3 and a degree 4 invariant. Note also that the numerator suggests a 
triquadratic syzygy in degree 6. The complete system of covariant polynomials was found by Le Paige in |31| . 
The simplest covariants is the ground form A. The three quadratic forms are 



and 



Bx[x) — det 
By{y) = det 

B^{z) = det 



d^A 

dxidzj 
d^A 



Q<i,j<l 



dxidyj 



/ 0<ij"<l 



To obtain the trilinear form, one computes anyone of the three Jacobians of A with one of the quadratic forms, 
which turn out to be the same 

dA dA 
dxQ dxi 
dB^ dB^ 
dxQ dxi 



C{x,y,z) = 
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The three quadratic forms B^^ By and B^ have the same discriminant A which is also the hyperdeterminant 
[Ti^ of the form. Furthermore the syzygy is 

+ ^B^ByB, + AA'' = 0. 

With each form A, we associate the vector va ■— {[Bx], [By], [B^], [C], [A]) where [P] = is P = and [P] = 1 
if P 7^ 0. The evaluation of va allows us to distinguish the different orbits (see Table [S]). 



Orbits 


Representatives 


VA 


Ovi 


|000) + 


(1,1,1,1,1) 


Ov 


|001) + |010) + |100) 


(1,1,1,1,0) 


Oiv 


+ |001) 


(0,0,1,0,0) 


Oiii 


+ |100) 


(1,0,0,0,0) 


Oh 


+ |010) 


(0,1,0,0,0) 


Oi 


1111) 


(0,0,0,0,0) 



Table 8: The case 2x2x2: evaluation of va on the orbits. 



Remark 5.1. Let us compare the orbits described in Tableland the covariants of Tabled Recall that in 
the 2x2x2 case, the variety of separable states, corresponding to O/ is X = x x P^. The invariant A is 
the equation of the dual variety X* , the so-called Cay ley hyperdeterminant but it is also the defining equation 
of the tangential variety t{X) according to Table [T] Its singular locus corresponds to hypermatrices I^P) such 
that C{\^))^Q. 

The syzygy restricted to hypermatrices |4') which belongs to t{X) becomes — —-B^ByB^- It tells us 

that the locus defined by C = is not irreducible but will be made of three components corresponding to the 
vanishing of one of the covariants B^, By, B^. To get a better understanding of the covariant B^ let us consider, 
for a given state l^f), the projective map -02; : P^ — >■ P^ defined by V'a, : — > (8) with 

r / X / ^0002^0 + ^1002^1 ^ooia^o + ^ioia;i\ ^ _ f 

V^oio2^o + ^1102^1 Aquxo + AuixiJ ~ J 

Let E = P^ X P^ c P'^ be the hypersurface defined by det = (the projectivization of the set of matrices 
of rank one in P(C^ ® C^)). The definition of the covariant B^ implies that Bx{\'i')) = if and only if 
ipxiV^) is contained in E. The image ■0a; (P^) is either a point or a line. If it is a point then 'tp^i'C'^) = f ® 9 and 
l^i) = [e®f®g\ e P^xP^xPi. If 0;,,(Pi) is a line then either 0;j;(pi) = P^ = ¥{C'^®g) or0^(pi) P^ = P(/®C2) 
(the variety S is ruled by two famihes of lines). The first solution implies I^E") G 0-2 (P^ x P^) x P^, i.e. \^) € O/y, 
the second solution gives I^E") S cr2(P^ xV}_x P^) x P^, i.e j^I') g O//. Therefore we recover geometrically that 
Bx{\^)) = <=> I*) G (ct2(F^ X F X pi) X pi) U (cr2(P^ X pi) X pi). Similarly it can be shown geometrically 
that By wiU vanish if and only if |*) G (P^ x cr2(P^ x P^)) U (cr2(P^ x P^) x P^) and B^ will vanish if and only 
if 1^') e (cr2(Fi X F X pi) X pi) U (pi x (T2(Fi X pi)). 

Remark 5.2. We also give an interpretation in terms of pencils of circles in appendix [BI 



5.3 The case 2x2x3 

From we find that the Hilbert series of the algebra of covariants is 



1 - 



1 - 



1 - 



xn. 



1-^ 

vfU3 



1 _ 

V3U3 
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The extraction of the positive part is a very tedious process which can be performed by expanding in simple 
fractions and erasing the fractions whose expansion as series are not Taylor. The Hilbert series of the algebra 
of the covariants can be explicitly computed : 



(1 - tWiU2M3)(l - t^uDil - t^UiU2){l - t4uiu2)(l - f^UiuDil - i^). 



(3) 



This suggests that there are 6 generators of degree 1, 2, 3, 4, 4 and 6 with only one invariant in degree 6 and a 
syzygy in degree 8. 

The covariant of degree 1 is just the ground form 

1 2 

^ = ^ ^ ^ ^ O'iJ.kXiUjZj-. 
i=0 j,k=0 

The covariant of degree 2 is obtained from A by elimination of variables x and y: 



B det 



dxidyi 



0<i,j<l 



The covariant of degree 3 is the bilinear form 



C 



aooo 
aooi 

xiyi 



OlOO 
Oioi 

ai02 
-xoyi 



aoio 
aoii 
aoi2 



am 
aii2 




0'i,j,cr{l)-l (ao,0,cr(2)-l,ai,l,cr(3)-l ~ '^l. , 0, a {2)- 10-0.1,(7 {3) -l) 1 



-xiyo xoyo 

The unique invariant generator is the determinant of C seen as a 2 x 2 matrix: 

A := det 



dx^dyj 



To describe the two covariants in degree 4, we recall the definition of the transvection of two multi-binary forms 
on the binary variables x^^^ = {x'^\ x^'i^ ),..., x^^'^ = {x'^\x'f^)): 



Ui9)ii,...,i 
where is the Cayley operator 



.f]X,/(:.'«,...,:z:'(P))g(x"W,...,x"(P)), 



d 
Ox" 

dx\ ^ 

J- OK, 







and tr sends each variables a;',x" on x (erases ' and "). 

The covariant A and C are two bilinear forms on the binary variables (xq, a;i) and (yo; 2/1 )• So we can apply the 
transvection operators and obtains two covariants in degree 4: 



Dx = (A, C)oi and By = (A, C) 



10- 



The evaluation of the vector va ■= {[B], [C], [D^], [Dy], [A]) on the different orbits is reproduced in Table [5] 
Note that the covariants and Dy have no role and that va has the same evaluation for each orbit in 
{Ovi, Ov, On} and for each orbit in {Ojv, Om, Oj}. So the knowledge of the covariant polynomials does not 
allow us to decide to which orbit a given state belongs. 
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Orbit 


VA 


Oviii 


(1,1,1,1,1) 


Ovii 


(1,1,1,1,0) 


Ovi 


(1,0,0,0,0) 


Ov 


(1,0,0,0,0) 


Oiv 


(0,0,0,0,0) 


Oni 


(0,0,0,0,0) 


On 


(1,0,0,0,0) 


Oi 


(0,0,0,0,0) 



Table 9: The case 2x2x3: Evaluation of va on the orbits. 



We need to compute more concomitant polynomials. For ternary variables, one has to consider a ternary 
contravariant variable C = (Co, Ci, C2) ^nd use an adapted version of the transvection: 

(/, 5, = ir^r,nln1^llf{x', y\ z', (^", y", C)Hx'\ y"\ z'", C") 

where 

d d 

"p — _a_ _a_ 

for p — X, OT p — y, 



a 


a 


a 


%" 


%° 


Of 










!li 


dv'2 


dp'2 


dp'2' 



a p = z or p = C and tr is the mapping which erases the symbol ', " and "', as previously. 
We define three concomitants: 

(A An)?,o,i' Byc ■■= (A A^c)o,i,i and := (S, Pc)%,o, 

where 

2 

i=0 

Let = {[B], [Bxc], [By(], [C], [A], [I?^]), we resume the evaluation of wa on the various orbits in Table [TUl 



Orbits 


WA 


Oviii 


(1,1,1,1,1,1) 


Ovii 


(1,1,1,1,0,1) 


Ovi 


(1,1,1,0,0,1) 


Ov 


(1,1,1,0,0,0) 


Oiv 


(0,1,0,0,0,0) 


Oiii 


(0,0,1,0,0,0) 


On 


(1,0,0,0,0,0) 


Oi 


(0,0,0,0,0,0) 



Table 10: The case 2x2x3: evaluation of wa on the orbits. 
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Remark 5.3. In the 2x2x3 case, the orbit Oj is the Segre product X = P-^ x P-^ x P^. The invariant A is the 
hyperdeterminant in the sense of of format 2x2x3, i.e. the equation of the dual variety X*. According 
to Table [21 the invariant A can also be interpreted as the equation of J{X,Oiv)- Like in the 2x2x2 case, 
the covariant C vanishes on the suborbits of the hypersuface defined by A = 0, i.e. C{\'^)) — means |\E') is a 
singular point of A = 0. The covariant C admits also the following interpretation in term of secant varieties: 

C{\^)) = ^ I*) e (T2(P^ X pi X p2) 

To prove this assertion, let us consider j^*) as a projective map -0 : P^ — > P'' defined by the linear map 

: ^ ® given by ^i{v) = Mq (g) e^{v) + Mi (g) el{v) + M2 (g) e^{v), with Mk = {Aijk)o<^,J<l e (g) 
for A: = 0, 1, 2 and e* the dual basis of (see remark [??T|l . The linear map ■(/; is of rank 3 precisely when C 7^ 0. 
But the rank of l^") S C"^ and the rank of -0 satisfy by construction 

rankd^f)) > rank(0) 

If I^P) is a limit of (the projectivization of) tensors of rank less than 2, by continuity of C we have Cd^*)) = 0, 
i.e. I^-) e (T2(P^ X Pi X P2) implies Cd^-)) = 0. 

On the other hand if Cd^*)) = we can assume without loss of generality that i/j = Mq (E) Bq + Mi (g) e*. If 
Mo = XMi we can write the map ?/; as -0 = Mq Bq, i.e. j^*) = (cq ® cq + ei (g) ei) (g) cq and |^) is a point of 
(72 (P"'^ X P-'^) X P^. If Mq and Mi are not colinear, these two matrices define a line after projectivization, i.e. 
pi = P(AMo + Aiil/i) C P''' ^ P(C2 (g) C^). Again consider S = pi x P^ the hypersurface of ¥{C'^ g) C^) defined 
by det = 0. If det(AMo + /iMi) = 0, i.e. the projectivized line defined by Mq and Mi is contained is S, then 
we can assume either Mq = bq g) bq and Mi = eg g) ei or Mq = cq g) cq and Mi = ei g) cq (both matrices are of 
rank 1 and their linear combination is of rank 1) and therefore either j^*) = eg g) cq (g eo + eo (g ei ® ei and \'^) 
belongs to P^ x o-2(P^ x P^) or l^f) = eo g) eg (g cq + ei g) cq (g ei and [*) belongs to cr2(P^ x F|_ x P^) x P^ (hke 
in example 15.11 the two possibilities correspond to the two families of lines of S = P^ x P^). 
If det(AA/i + 11M2) ^ for some values A, /i, i.e. the line intersects E, then we can assume Mq is a 2 x 2 matrix 
of rank 2 and Mi a matrix of rank 1. There are two cases to consider: 

• The line is tangent to S and we can assume Mq = eog)ei+eig>eo and Mi — eo g> cq. Then |^) = 

(eo g) ei + ei g) eo) g) gq + eo ® eo g) ei = eo g) ei g) eo + ei g) eo g) eo + eo g) eo g) ei = ^(— eo g) eo g) eo + 
(eo + eei) g) (eo + eei) g) (eo + eei)). From \^)^ = ^(-eo g) eo g) eo + (eo + eei) g) (eo + eei) g) (eo + eei)) 
we see that is a limit of tensor of rank 2 when e — > 0, i.e. is in the closure of cr2(P"'^ x P-'^ x P^). 

• The line is secant to S and we can assume Mq = eog)eo + eig)ei and M2 — eo g) eo. Then \^) = 
eo g) eo g) eo + ei g) ei g) ei, i.e. I^*) S cr2(F^ x P^ x P^). 

The covariant B is the analogue of the covariants Bx,By,Bz in the 2x2x2 case, i.e. B{\'i>)) — det(0(C^)). 
The vanishing of B{\'i>)) implies ?/'(P^) belongs to E. In particular one sees that B vanishes only if C = (C 7^ 
implies '0(P^) is a plane and therefore can not be contained in E). Like in the 2x2x2 case there will be three 
cases corresponding to i? = which are: 

• ?A(p2) is a point of E (orbit O/ = P^ x P^ x P^), 

• ^(F^) is a hne of E, i.e. 

- either 0(P^) = P(/ ® C^) (orbit Ojv = P^ x (73 (P^ x P^)) 

- or ?A(P2) = P(C2 g) g) (orbit Om = cr2(P^ x P^ x P^) x P^). 

Those different cases are distinguished by the concomitants B.j.,q and ByQ. If C = and B ^ 0, then there are 
three different cases 
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• i/>(p2) is a point and V(P^) i S (orbit On) 

• 7/'(P^) is a line tangent to S (orbit Oy) 

• '(/'(P^) is a line secant to S (orbit Ovi) 

Finaly when C ^ then ?/;(P^) is a plane which could be tangent to E, A = (orbit Oy//) or secant, A 7^ 
(orbit Oviu). 

Remark 5.4. In appendix |B] we give an interpretation in terms of linear complex of circles. 
5.4 The case 2x3x3 

Again we start by computing the multivariate Hilbert series: 

r!"^CT„,l]t^CT„,r!"^ ( 1 _ 1 Vl - ( 1 - ( 1 - ( 1 - (l-- 




^ ^ ^ ^ -^'^ V2ul)X V3UIJ\^ U2vV\ ^3^1/ \ ^2 

Ul H I \U2+V2-\ ) \U3+V3 + 



x((i-^)n, 

U3 



Ul J \ V2U2 I \ V3U3 



Wc find 

P{t;ui,U2,U3) 



(4) 



(1 - tUiU2U3){l - Ult^){l - U2U3t^)il - ult^){l - _ ulult^){l - u2f6)((l _ t^2) 

with 

P{t; Ul, U2, U3) = -t26u^4y25yg5 _ ^22y^2y24y^4 _ ^Zly^Sy^Sy^S _ il9u^3y24y^4 _ t^Sui^U2^U3^ 

-ii Vu2^Mi - (mi - 1) (ui + 1) t^^ua'^ua^ - (ui - 1) (ui + 1) t^^U3^U2^ui 

- (Ui - 1) (Ui + 1) i"u32u2^Ui - (mi - 1) {Ui + 1) tl°M3U2Ul^ 
+ t^U3^U2^ + fu3U2Ui + t^W3^M2^Ui + t*U3U2Ui'^ + 1. 

This suggests a very complicated description of the algebra. Nevertheless, as in the other cases, we will use 
only a part of the concomitant polynomials. The degree 3 generator is a cubic binary form obtained from the 
ground form A := X^Lo Y.].k=Q°-i3kXiyjZk by: 

Cx ■= det 



For a cubic binary form, the algebra of covariants is well known. Its Hilbert series is 



1 - uft'^ 



(1 - tuf){l - ujt^){l - ult^){l - t^) ■ 
Considering a generic cubic binary form 

a := ooXq + aixixl + a2xlxo + asx^, 
we compute a covariant of degree 2, the Hessian: 

b := det f jj^) 

= (3aoa2 - a\)xl + (9aoa3 - aia2)xiXo + (80103 - 02)^X1. 

The discriminant of this quadratic form is the only invariant generator: 

d := 4(3aoa2 — a^)(3aia3 — a\) — (9aoa3 — 0102)^. (6) 



(5) 
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We need also to compute the covariant of degree 3: 

c := (6 zo^o ^2 — 2 zoflj ^ — ij -Zifla + 9 flo 03 zi) (aj zq^ + 2 flg ziZq + 3 03 zi^) 

— (-zoflj a2 +9 zqUo 03+6 ziBj 03 — 2 zia^^) (3 co zo^ + 2 zizo + aa zi^) . 

Note that we have the syzygy: 96 + 512c + 128^0^ = 0. 

Replacing each at by the coefficient of XqX^^' in Cx in eq. (|5l6l7p . that is, 











oooi 


O002 












ao 






OOll 


O012 














ao20 


0021 


O022 










flooo 




aoo2 






oooi 


O002 




oioo 


OlOl 


O102 


ai = 


floio 


OOll 


aoi2 


+ 




Olll 


0112 


+ 


oqio 


Ooll 


O012 




0-120 


ai2i 


ai22 




1020 


0021 


O022 




0020 


O021 


O022 




flooo 


flooi 


aoo2 




OlOO 


OlOl 


O102 




Oloo 


OlOl 


O102 


02 = 


fliio 


am 


aii2 


+ 


fliio 


Olll 


0112 


+ 


oolo 


ooli 


O012 




0120 


ai2i 


ai22 




0020 


O021 


O022 




O120 


0121 


0122 








flioo 


OlOl 


O102 
















Alio 


Olll 


0112 














ai20 


0121 


0122 









we obtain one quadratic covariant in degree 6, a cubic covariant in degree 9 and the invariant A in degree 
12. 

We observe that = if and only if I^; = 0, so we have only to evaluate the vector va — {[Cx], [Fx], [A]} on 
the representative of the different orbits. These covariants are not sufficient to discriminate between the orbits 
but they allow us to partition the set of the orbits into 4 sets (see Table [TT|) . The other covariant polynomials 



Orbits 


VA 


Oxvii 
Oxvi, Oxiv 
Oxv, Ox III, Oix 
Oi, . . . , Oviii, Ox,Oxi, Oxii 


(1,1,1) 

(1,1,0) 
(1,0,0) 
(0,0,0) 



Table 11: Evaluation of va on the orbits: the case 2x3x3. 



are more complicated to understand in an algebraic way. For instance, we will use the two following degree 4 
covariants: 

Dxyz {Cx, A) 

and 

Dy, := tr {nxtv'xn',n,nyA{x',y', z'fA{x", y" , z")A{x"\ y'" , z'")) 



where $7' :— 



a d 



and trf, sends x'- and xT to x'- . The covariants Dyz and D^,,z are both bilinear in the 



'xyz 



Fix" " 

ternary variables y = {y^, yi, j/2} and z — {zq, zi, Z2} and Dxyz is quadratic in the binary variable x — {xq, xi}. 
These covariants are used to separate the orbits Oxv, Oxiii and Oix (see Table [T^). Define also 



trQx^^Aix', y', z')A{x" , y" , z")Dy,{y"' , z'"). 
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Orbits 




Oxv 
Oxiii 
Oix 


(1,1) 

(0,1) 
(0,0) 



Table 12: Evaluation of {[Dyz], [Dxyz]) on the orbits Oxv, Oxiii and Oix- 

This polynomial vanishes on Oxiv but not on Oxvi- It remains to separate the orbits O/, . . . , Oviii, Ox, 
Oxi andOxii- 

We need to introduce some covariant polynomials: 



:= trnyn^Aix', y\ z')A{x" , y", z")P{y"' , ry'"), 



and 



By^ := trfl,f7,A(x', y', z')A{x\ y" , z")P{z'" , ('")■ 
These polynomials allows to discriminate between the orbits Oj, On and Om (see Table [T51) . In degree 2, we 



Orbits 


{[BycliBz^]) 


Oin 
On 
Oi 


(0,1) 
(1,0) 
(0,0) 



Table 13: Evaluation of {[By^], [B^r;]) on the orbits Oi, On and Om. 
define also the concomitant: 



In degree 4: 
and 

Also in degree 6: 

Finally, in degree 8: 
and 



B,,C := iA,A,P{y,r^)P{zX)r'. 
Dyzr,Q tr^ynz(n,,A{x' , y', z')A{x" , y\ z")fP{y'\ v"W, C'" )■ 

{ByQ, ByQ, Bycf^"^ . 

Hxyzv ■■= trn^Dyz^iiy', z' , 77', C)Bz^{z", ri")Bx^^{x"' , 77"', C'") 



Hxy.C tvnr^DyzrKiy', z\ vj , OBydv", C) B,,d^'" , v'" , C") 

We define v'^ := {[Dyz-rjc], [Bxr^i^], [Dnc], [Dy^], [Fr^, [Fq], [ff^azcl, [HxyziT\)- TablelHgives the evaluation of on 
the remaining orbits. 

We summarize the results of this section in Table (TSl setting 



WA {[Dyzr^c], [Bx-nc], [By(^], [B^^], [Cx], [A,c], [Dyz], [Fy], [Fx], [F,^], [Fq], [Hxyzc], [Hxyz-q], [A]). 
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Orbits 


v'a 


Oxii 


(1,1,1,0,1,0,0,1) 


Oxi 


(1,1,1,1,0,0,1,1) 


_Ox 


(1,1,1,0,0,0,0,1) 


Oviii 


(1,1,1,0,0,1,1,0) 


Ovii 


(1,1,1,0,0,0,1,0) 


Ovi 


(1,1,1,0,0,0,0,0) 


Ov 


(1,1,0,0,0,0,0,0) 


Oiv 


(0,1,0,0,0,0,0,0) 



Table 14: Evaluation of on the orbits C/y, . . . , Oviii, Ox, Oxi and Oxii- 



Orbits 


WA 


Oxvii 


(1 


1, 


1 


1 


1 


1 


1 


1, 


1 


1, 


1 


1 


1, 


1) 


Oxvi 


(1 


1, 


1 


1 


1 


1 


1 


1, 


1 


1, 


1 


1 


1, 


0) 


Oxv 


(1 


1, 


1 


1 


1 


1 


1 


1, 





1, 


1 


1 


1, 


0) 


Oxiv 


(1 


1, 


1 


1 


1 


1 


1 


0, 


1 


0, 





1 


1, 


0) 


Oxiii 


(1 


1, 


1 


1 


1 


1 





0, 





0, 





1 


1, 


0) 


Oxii 


(1 


1, 


1 


1 





1 





0, 





1, 








1, 


0) 


Oxi 


(1 


1, 


1 


1 





1 


1 


1, 





0, 





1 


1, 


0) 


Ox 


(1 


1, 


1 


1 





1 





0, 





0, 








1, 


0) 


Oix 


(0 


1, 








1 








0, 





0, 








0, 


0) 


Oviii 


(1 


1, 


1 


1 





1 





0, 





0, 


1 


1 


0, 


0) 


Ovii 


(1 


1, 


1 


1 





1 





0, 





0, 





1 


0, 


0) 


Ovi 


(1 


1, 


1 


1 





1 





0, 





0, 








0, 


0) 


Ov 


(1 


1, 


1 


1 











0, 





0, 








0, 


0) 


Oiv 


(0 


1, 

















0, 





0, 








0, 


0) 


Oiii 


(0 


0, 


1 














0, 





0, 








0, 


0) 


On 


(0 


0, 





1 











0, 





0, 








0, 


0) 


Oi 


(0 


0, 

















0, 





0, 








0, 


0) 



Table 15: Case 2x3x3: Evaluation of wa on the orbits. 

Remark 5.5. The comparaison between varieties of Table HI and [T5] give the following interpretations for A, 

Again the only invariant polynomial A can be considered as the equation of the dual of X = x x or 
as the equation of the join J{X, t{X)). 

The covariant is defined by Ca;(|*)) = dct(?Aa;(C^)) where i/j^ : ^ ® C^, is given hy ip^ — (E) Mq + 
el ® Ml, with Mi 3 X 3 matrices. Thus C^H^i')) ^ means tpxiP'^) is not contained in E = cr(p2 x P^) c P^, 
the hypersurface defined by det = 0. This is the case when j'i') is a general point of Ojx — x P®. Thus Cx 
does not vanish on all orbits containing Ojx, i-e. orbits Oxiii,Oxiv,Oxv,Oxvi,Oxvii (see FigureH]). On 
the other hand when Cxi\1')) = 0, then Vx(F^) is a subset of S = dailP^ x P^) c P^ 

The image ipxO^^) is either a point or a line. It is a point when By(^{\'^)) = and Bzr^H'^)) = i.e. the matrices 
Mo and Mi which define tpx are colinear. The discussion on Cx and Sy^, B^ri gives the following possibilities: 

• if By^il^)) = B.rjiM) = iM^') is a point of P«) 

- then C:r(|*)) = (i.e. ^(IP^) e S), gives two possibilites: 
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* is a point of x (orbit Oi) 

* Va;(IP^) is a point of S (orbit Oiv) 

- and Cxi\^))^Q (i.e. i/'(P^) ^ S), gives one possibility: 

* Va;(]P^) is a general point of \ E (orbit O/x). 
• if Byc^{\'^)) ^ or B^,,(|*)) ^ (V'x(P^) is line in P**) 

- then C^{\^)) = (i.e. ■(/'(P^) C S), gives 9 possibilites: 

* Vx(pi) is a line of P^ x (orbits O// and 'Oui) 

* Vx(P^) is a line tangent to P^ x P^ (orbit Oy) 

* is a line secant to P^ x P^ (orbit Ovi) 

* Va;(P^) is a line of E (orbits Oyiii, Ox, Oxi, Oxii) 

- and C^d^-)) 7^ (i.e. -(/-(P^) n E = 0), gives 5 possibilities: 

* V.(PI) is a line of P* \ E (orbit Oxiii, Oxiv, Oxv, Oxvi, Oxvii)- 

The covariants and concomitants allow us to distinguish the different positions of the lines but a priori geometric 
interpretations of those polynomials are far from being obvious. 

Note Table [TSl have been computed using Maple programs. 

The sources are available at |http : //www-igm. univ-mlv . f r/^luque/f orm233 . txt 

6 Conclusion 

In this paper we proposed an alternative approch to the geometric descriptions of entanglement given by 
Miyake in |37l 139) . The idea was to use auxiliary varieties such as join, tangent and secant varieties instead of 
a description by dual varieties. The introduction of the secant and tangential varieties brought a more precise 
description of the singular locus of the dual varieties as we were able to interpret the singular components 
of X* as dual varieties of the stratification by join and tangential varieties. We also detailed the geometric 
description of the entangled states for the 2x3x3 quantum system. Both descriptions of entanglement, 
by join and tangential varieties or dual varieties, are equivalent as long as we deal with group actions with 
finetly many orbits. However challenging problems in QIT start with quantum systems with infinitely many 
orbits. For example in the case of the 3x3x3 system or the 2x2x2x2 system there does not seem 
to be a complete consensus on what an entangled state is mathematicaly (different papers announce different 
numbers of entangled states under different definitions |49|[T1[27]). We believe that the approach by secant and 
tangential varieties could bring interesting perspectives for the geometric description of entanglement in QIT. 
For instance the recent work of Buczyhsky and Landsberg on the third secant varieties of the Segre product of 
three projective spaces provides useful results to describe the points lying in the closure of 173 (P'"'! x P'^^ x f''^^). 
But smooth points of cr3(P''i x P'''^ x P''^) correspond to a state of type |000) + |111) + |222), i.e. generalize the 
GHZ state. As noticed in [16j the results of [S] on the geometry of (Ts(P^ x ■ . ■ x P^) should lead to a better 

n times 

understanding of the quantification of entanglement for n-qubits. 

Another way to study mathematicaly entanglement is to look for invariant polynomials under the SLOCC group 
actions [H [M] . These polynomials allow one to distinguish the different states as we did in Section [S] where we 
introduced an algorithmic method to identify the orbit of a given state for the quantum systems of this paper. 
In this sense the dual varieties approach provides also invariant polynomials as the hyperdeterminant, in the 
sense [H], is SLOCC invariant. But the study of the defining equations of the secant and tangential varieties is 
also a topic of interest for algebraic geometers [HI |^ . Looking for geometric interpretations of the covariants we 
made some connection between the vanishing of some these polynomials and the defining equations of the secant 
variety. We believe that the geometric understanding of the polynomials produced by classical invariant theory 
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techniques in the context of auxiUary varieties needs to be better understood. The study of the geometry of 
secant and tangent varieties of Segre varieties should therefore help to understand the structure of multipartite 
entanglement both from the orbit stratifications point of view and from the invariant polynomials perspective. 



A Hilbert series 

In this section, we use extensively symmetric functions (and in particular Schur functions) and their applications 
to the representation theory of the linear groups (see e.Q.|35|). 

The characters of the irreducible polynomial representations of the product group G are product of Schur 
functions 

where the A*-*-* are partitions and each Xi = {xn, . . . , Xim} is a set of Ui variables. The character of the 
one-dimensional representation 

det{gY ^det{g,Y' ...detigkY". 

is the product of rectangular Schur functions s^^i (Xi) . . . s^^k (^fc), whilst the character of G is Sd{Xi . . . Xk). 
Hence, the dimension of the space of invariants of degree d and weight i, which is also the multiplicity of the 
one dimensional character dct^ in S'^{'H), is given by the scalar product 

dimlnv{d, i) = {sd{Xi . . . Xk), s^i^i (Xi) . . . s^k iXk))G 

of characters of G. To evaluate this scalar product, we can replace the Xi by an infinite set of independent 
variables, and compute in the tensor product Sym®^ of k copies of the algebra of symmetric functions Sym. The 
results will be the same, since in both cases the orthonormal basis is given by tensor products of Schur functions, 
S\ being identified with Sx^i) ® ■ ■ ■ ® Sxw ■ Under this identification, the operation 5{f) = f{XY) corresponds 
to a comultiplication in Sym, which is known to be the adjoint of the internal product * of symmetric functions. 
Note the value of £i depends on whose of d and rii {iiUi = di). Hence, 

dimlnw(d) = {sd, * • • • * )sym (8) 

A similar reasoning gives 

dimCou(do, ■■.,dk) = {sdo, (sri^di) * ' • ' * Sdfc))sym, (9) 

1 k 

again the value of each £i is obtained by di + liUi = do- 



In order to compute the Hilbert series we need to introduce the Cauchy function Ht which is a very powerful 
tool for the manipulation of symmetric functions (see e.g. [29j): 



xeX [n>l J ri>0 



where p„ — X^^jgx ^" denotes a power sum symmetric function. 
Consider also the operator 

4:=expJ-^(-t)" ^ 



dPn 



and the vertex operator |30| : 

Tt ■.= Tltd_i, 
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The operator dt acts by shifting the power sum 

dtf{pi,P2, . . . ) - fiPl +t,P2-t^,...,). (10) 

The definition of Schur function can be naturally extended to the compositions with negative parts [29) 



5^2 + 1 



5t?i+n — 1 



(11) 



with s„ = is n < 0. Using equahties pU|) and (|lip we find 
And we use it iteratively to show 



(12) 



where YX = {yiX : i — I . . .n^x G X}. 

On the other hand the dimension formula ^ allows us to write the Hilbert series as a scalar product 
involving Cauchy functions. Let us denote by [] the plethysm operation (see e.g. |35j). Remarking that for an 
alphabet of size k we have: 

SlkSd = ^ S„[si + Sifc], 
l,d n 

we obtain 

Hcov{t;u) = dim Gov {do,..., dk)t'^"u^' ...u^" = (Hi [a.^ (t, ui)], Hi [a,^ (1, iia)] . . . ni[ai, (1, 

(13) 

where ak{t, u) — tusi + t'^Sik for fc-ary alphabet. Note we have 

A = (Ai,...,Afc) 



this is a consequence of 



A = (Ai....,Ap) 
Ai + ... + Afc = „ 



and of the Fieri formula. Hence, for an alphabet of size k we have 



>-=(>-l >-k) 



These series can be obtained using a combination of vertex operators and Omega operators of Macmahon 
ri> [36J which send the monomials with a negative power of u to 0. Indeed, from ((T^) we obtained for an 
alphabet of size 2 (binary case) 



e{\)<2 
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and for an alphabet of size 3 (ternary case), by setting yi — u, y2 — and = v in (|12p 

- \ VW^ ) V U) \ V^U J 

where CT^, means the constant term. 

Combining this with equahty p3p this gives the Hilbert series as an Omega 

Hcov{t- u) = CT,„,r!^i . . . CT,„ri^"B,, {ui,vi)... (u„, z;„)nt [A,, (ui, vi) • • • w„)] (14) 

where ^2(1*, w) = u + i, i?2(w, u) = 1 — for binary variables and 

A2(u,w) =ii + w + — , B3(m,v) = (1-^) (1--) 

for ternary variables. 



u^ h I X 

vu 



B Interpretations in terms of annalagmatic geometry 

Let us recall first the principle of tetracyclic coordinates [TT] . 

One starts with the equation of a circle in the real Euclidean plane with coordinates {x,y). Introducing 
homogenous coordinates {X :Y:Z),on the complexified plane, the equation of our circle C can be written in 
the form 

xo ■ i{X^ + Z^) + xi ■ (X^ - Z^) + X2 ■ 2XZ + X3 ■ 2YZ = , (15) 

where {xq : xi : X2 ■ X3) are called the (homogeneous) coordinates of C. The quantities yo = i{X'^ + Y'^ + Z^), 
2/1 = X^ + Y^ — Z^ , y2 — 2XZ, y^ — 2YZ, are called the special tetracyclic coordinates of the point (X : Y : Z) 
of P2. They satisfy 

(yy) = o, (16) 

where 

3 

(xy) = ^ x^yi (17) 

i=0 

is the fundamental quadratic form in the geometry of circles. Any other nondegenerate quadratic form can be 
taken instead of (xy). We shall set for a circle 

Xq = a + d , xi = i{a — d) , X2 = b — c X3 — i{b + c) , (18) 

so that 

(xx) = 2(ad- 6c) = 2det M {xx') ^ ad' + ad ~ be - b'c (19) 

where M is the matrix 

M=(:2). (20) 

The tetracyclic coordinates of a point is now given by a rank 1 matrix which can be written in the form 

a ^\ ^ ( 2/2^2 -y2Zi 
c ~ \~y1Z2 yizi 



(21) 
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In this picture, finite points of the plane are represented by points on the quadric S: detM — 0, and the 
coordinates of a point are given in terms of the parameters (^1,2/2) and (21,22) of the two generatrices of the 
quadric intersecting in it. The finite points are those for which ixo + xi = —2Z^ ^ 0. In matrix coordinates, 
this reads a ^ 0. Hence, the elements at infinity of this geometry consist of the two intersecting generatrices 
b = and c = of S. We are working with a compactification of the complexified Euclidean plane which 
is not isomorphic to but to x P^. The generatrices of E correspond to two families of lines of the affine 
plane, called minimal lines. 

Turning back to the multilinear form 2x2x2, we can now interpret the equation X)i j=o ^ij^iVj = as a that 
of a circle, whose coordinates form the matrix (a^ ). Hence, the equation A = can be understood a describing 
a pencil of circles. To do this, we have to single out one variable, say z, viewed as a projective parameter, the 
other two ones being minimal line coordinates on a tetracyclic plane. In this setting, the orbit classification and 
the normal forms are almost immediate. Let Cz be the point in P'^ representing the circle of the pencil with 
parameter z. Then, the set £ — {Cz\z £ can be either a proper line (the generic case) or be degenerated into 
a single point. The rest of the discussion will depend on the relative position of this line or point with respect 
to the non-singular quadric E of null circles. Let us first consider the case where £ is a proper line. Generically, 
it will intersect E into two distinct points Ci, C2 (the base points of the pencil). These points can be mapped 
by a circular transformation to the origin of the affine plane, and to the intersection of the two isotropic lines 
at infinity, which have respectively as matrix coordinates 

' and (» »■) (22) 



oy \o 1^ 

whence the normal form associated to Oyi in this case. The test for this case is A ^ 0, since = gives the 
parameters for the two null circles of the pencil. We note on the normal forms that C = is a circle of the same 
pencil. It is the unique such circle which is orthogonal to Cz- Let us remark that since we are working over C, 
we do not distinguish between intersecting pencils and pencils with limit points, but it would be easy to refine 
the discussion in order to obtain the classification over R, which can be found in Sokolov's book |44| . 
The next case is when i is tangent to E in exactly one point. Here, A = 0, and there is exactly one null circle 
Cz- By a circular transformation, we can arrange that this null circle becomes the origin, and that the radical 
axis becomes a coordinate axis, say the y axis whose equation in matrix tetracyclic coordinates is 6 — c = 0, 
d = 0, whence the normal form associated to Oy for this case. All the circles of the pencil have its null circle 
as a common point. Its equation is given by C = 0. 

By a further degeneracy, £ can become a generatrix of E. The two systems constitute then two orbits, with 
respective normal forms associated to On and O///. 

Finally, £ can be a single point m. If it is not on E, we can transform it into any proper circle, e.g., the one 
with matrix coordinates 

^0 1' 

which yields the normal form associated to Ojv- If m is on E, it is a null circle, and the normal form can 
obviously be taken as whose associated to Oj. 

For the case 2x2x3: the equation A = represents then a linear complex (or net) of circles Cz, i.e. a 
plane of the P^ of circles, which may degenerate into a line or a point. Let, as above, H^, be the linear variety 
formed by the representative points of all these circles. The ground form, identified with the matrix Mz = 
(o-ijozo + ciijizi + aij2Zo) - ^ = MiZi, is then the matrix of circle coordinates of C^, and B = det(M2) = 
when Cz is a null circle. 

In the generic case. Mi, M2, are linearly independent, and we have a proper net. This case is recognized 
from the covariant C, which does not vanish identically. The net is then formed by the collection of circles 
orthogonal to a fixed circle, whose equation is C = 0. By a circular transformation, this circle can be mapped to 
a coordinate axis, say Oz, whence a simple normal form (the net of circles centered on Oz). This is not anymore 
possible if detC = 0. The net is the formed by all circle having a common point, which me may take as the 
origin of coordinates. 
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If C vanishes identically, the net reduces to a pencil or to a point, and the normal forms can be inferred from 
the previous discussion of trilinear forms. 

To conclude, let us remark that D^Dy = gives the minimal lines through the common points of ^ = and 
C = 0. Also, the hypcrdctrminant A is proportional to the discriminant of B, whose detorminantal expression 
is recognized as the condition that the three circles with coordinates Mj have a common point. 
Acknowledgments This paper is partially supported by the ANR project PhysComb, ANR-08-BLAN-0243-04. 
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